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Generalized Parton Distributions and 
the Spin Structure of the Nucleon 

Xiangdong Ji®- * 
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Generalized parton distributions are a type of hadronic observables which has recently stimulated great interest 
among theorists and experimentalists alike. Introduced to delineate the spin structure of the nucleon, the orbital 
angular momentum of quarks in particular, the new distributions contain vast information about the internal 
structure of the nucleon, with the usual electromagnetic form factors and Feynman parton distributions as their 
special limits. While new perturbative QCD processes, such as deeply virtual Compton scattering and exclusive 
meson production, have been found to measure the distributions directly in experiments, lattice QCD offers a 
great promise to provide the first-principle calculations of these interesting observables. 


Since the EMC experiment on polarized deep- 
inelastic scattering (DIS) on the proton in the 
late 1980s ||], QCD spin physics has evolved into 
one of the most active frontiers in hadron physics. 
The International Spin Physics Conference is held 
every other year with an average of 350 physi¬ 
cists participating in the meeting. The topics of 
the conference are often dominated by QCD spin 
physics. There are many experimental facilities 
around the world which have been playing a key 
role in spin physics research. For example, at 
CERN, following the EMC experiment there was 
SMC, and now the COMPASS experiment is on 
its way to take data |^. At SLAC, a series of ex¬ 
periments were finished in the 1990s: E142, E143, 
E154, and E155. At DESY in Germany, the HER¬ 
MES collaboration has taken data for a number 
of years, and new runnings have been planned 
(HERMES H). A summary of the finished ex¬ 
periments can be found in The spin physics 
at RHIC has just started, where 250 GeV polar¬ 
ized proton beams can be used to explore spin- 
dependent high-energy collisions |Q|. Jefferson 
Lab in Virginia has a high intensity 6 GeV elec¬ 
tron beam, and many experiments there involve 
polarized beams and targets. A plan has been 
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made to upgrade the beam energy to 12 GeV. 
MIT-Bates delivered a polarized electron beam a 
few years ago, and a number of intereting low- 
energy experiments have been done which have 
an impact on spin physics. 

There have been many important theoretical 
developments in spin physics as well. Perturba¬ 
tive QCD analysis of gi{x) and g 2 {x) structure 
functions of the nucleon has been carried out 
to two- and one-loop orders, respectively 
The role of the axial anomaly in understanding 
the EMC results has been explored and debated 
1^. A classification of the leading and higher- 
twist polarized quark and gluon distributions has 
been made [^. QCD angular momentum and its 
role in high-energy scattering have been clarified 
1^. The effects of parton transverse momentum 
in high-energy processes have been explored sys¬ 
tematically 0 - Generalized parton distributions 
and deeply virtual Gompton scattering are the 
main subject of this talk. 

One of the most important goals in QCD spin 
physics is to understand the spin structure of the 
nucleon, i.e., how the nucleon spin is made from 
its fundamental constituents. Because nonpertur- 
bative QCD tools, such as lattice QCD, is still un¬ 
der development, we have relied, for many years, 
on model descriptions. In the naive SU(6) quark 
model, the three valence quarks move in the s- 
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orbit, and the spin of the nucleon comes entirely 
from the coupling of the quark spins. There are 
a number of experimental evidences which seem 
to support this simple picture. For instance, the 
magnetic moment of the nucleon can be calcu¬ 
lated in terms of the SU(6) wave function, and 
the ratio between the neutron and proton is fill] 

/in/Mp = -2/3 . (1) 

The experimental data is —0.68, very close! 
Moreover, the A ^ TVy transition can take place 
with both Ml and E2 electromagnetic multipoles. 
The quark model predicts that the E2 to Ml 
ratio is 0. Experimentally, the result is 1-3%, 
depending on how the background contribution 
is modelled. [However, recent results from the 
large QCD seem to indicate that these quark 
model predictions follow simply from the con¬ 
tracted SU(4) symmetry present in the large Nc 
limit, independent of the detailed dynamics [p^ ]. 

This simple picture of the spin structure of 
the nucleon has been challenged by the polarized 
DIS data from CERN, SLAG and DESY §. In 
these experiments, the polarized quark distribu¬ 
tions can be extracted: 

Aq{x) = q+{x) - q-{x) , (2) 

where ^±( 2 :) is the density of quarks with lon¬ 
gitudinal momentum fraction x and spin aligned 
(anti-aligned) to the spin of the nucleon. More 
precisely, the measured gi{x) structure function 
of the nucleon is related to the helicity distribu¬ 
tions by 

9iix) = ^ ( 3 ) 

/ 

where / sums over light quark flavors. A sum¬ 
mary of world data on the gi(x) structure func¬ 
tions for proton, neutron and deuteron is shown 
in Fig. 1. 

The operator production expansion (OPE) 
yields the deep-inelastic sum rule 

J g{{x)dx = ^{4Au + Ad + As) , (4) 

where Aq is the fraction of the nucleon spin car¬ 
ried by the quark spin of flavor q. The total quark 
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Figure 1. A compilation of world data on the 
gi{x) structure function of the proton, neutron 
(^He), and deuteron. 


spin contribution to the proton spin 
f dx Aq{x) = AS = Au + Ad + As , (5) 

cannot be extracted from the DIS data alone. 
One must use the axial couplings gA determined 
from neutron and hyperon beta decay as well as 
the SU(3) flavor symmetry. An analysis of the 
world data to the next-to-leading order in per¬ 
turbative QCD yields [^ 

AS(1 GeV^) =0.2±0.1 . (6) 

Therefore, at least 70% of the proton spin resides 
in either the orbital motion of the quarks or glu¬ 
ons. 

To understand the spin structure of the nucleon 
in the framework of QCD, one must start from the 
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QCD angular momentum operator in its gauge- 
invariant form Ip^: 


'^QCD — Jq Jg 1 

where 


Jq = 


(fx X xTa 


2-1 —* 

ijj''' —tjj + 'ip^x X 


•h — 


(fix x X (E X B). 


( 7 ) 


( 8 ) 


The quark and gluon parts of the angular mo¬ 
mentum are generated from the quark and gluon 
momentum densities Tq and E x B, respectively. 
S is the Dirac spin-matrix and the correspond¬ 
ing term is the quark spin contribution which we 
have discussed in the context of the polarized DIS. 
D = V — igA is the covariant derivative and the 
associated term is the gauge-invariant quark or¬ 
bital contribution. 

Using the above operator, one can easily con¬ 
struct a decomposition for the spin of the nucleon. 
Consider a nucleon moving in the z-direction, and 
in the helicity eigenstate A = 1/2. The total he- 
licity can be calculated as an expectation value of 
Jz in the nucleon state, 

-^-/AE{p) + Lq{g)->rJq{g), (9) 

where the three terms denote the matrix elements 
of three parts of the angular momentum operator 
in Eq. (||) . The physical significance of each term 
is obvious, modulo the scale and scheme depen¬ 
dence indicated by /i. The scale dependence in 
is generated from the U(l) axial anomaly 
j^. Note that the individual term in the above 
equation is independent of the momentum of the 
nucleon. In particular, it applies when the nu¬ 
cleon is travelling with the speed of light (the in¬ 
finite momentum frame) |p^ . 

The DIS data AE « I indicates that the nu¬ 
cleon is much more complicated than the naive 
quark model depicts. At a phenomenological 
level, one may regard the quark model as an ef¬ 
fective description but it is of little use to explain 
the DIS data if one does not know the relation 


between the QCD and constituent quarks. The 
challenge we have now is to find additional con¬ 
tributions to the nucleon spin in the context of 
the above decomposition. 

From the definition, we see that Jq g are just 
the matrix elements of the moment of the energy- 
momentum tensor in a nucleon helicity state 





/ 


d^x{x X Tq^gY 



( 10 ) 


I will argue below that the matrix elements can 
be extracted from the form factors of the quark 
and gluon parts of the QCD energy-momentum 
tensor Tj^Yg in the nucleon state. 

To motivate the linkage, let’s consider an anal¬ 
ogous physics observable. In the electromag¬ 
netism, the magnetic moment is defined as a spa¬ 
tial moment of the current density 



( 11 ) 


If one knows the Dirac and Pauli form factors of 
the electromagnetic current, Fi{QY and F 2 [QY^ 


{P'\f\P) = U{P')[Fi{t)r 

+F2{t)' 


2M 


U{P) , (12) 


where t = [P — P'Y, and A = P — P'. The spin- 
flip form factor AGmY) = + F 2 {t)) yields 

the electric current distribution in the nucleon. 
From this, we can calculate the magnetic moment 


Mp = Gm (0) = Fi(0) + F2(0) , (13) 


in units of eh/2Mc. 

Using Lorentz and discrete symmetries, it is 
easy to show that the symmetric and traceless 
part of the quark and gluon energy momentum 
tensors each supports three form factors, 


{P'\TYYg\P) = U{P')[AqJth^>^P'''^ 

+Bq,g{t)P^^ia^'>‘^Ac./2M 

+Cq,g (<) A(^ A") /m] U{P) , (14) 


where P = (P -1- P')/2. Taking the forward limit 
of the /i = 0 component and integrating over 3- 
space, one finds that Aq^g{0) give the momen¬ 
tum fractions of the nucleon carried by quarks 
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and gluons (the momentum sum rule contraints 
Aq{0) + Ag(0) = 1). On the other hand, substi¬ 
tuting the above into the nucleon matrix element 
in Eq. ([l^), one finds ||^ 

“ 2 )] • (15) 

Therefore, the matrix elements of the energy- 
momentum tensor yield the fractions of the nu¬ 
cleon spin carried by quarks and gluons. 

While a measurment of the gravitational form 
factors is impossible in the forseeable future, they 
can be extrated from QCD hard scattering. A 
process which was first identified is deeply vir¬ 
tual Compton scattering | p^ , in which a space¬ 
like virtual photon created through lepton scat¬ 
tering, scatters off the nucleon target, yielding a 
final state of a real photon and a recoiled pro¬ 
ton (see Fig. 2). When the virtual photon is in 
the Bjorken limit, the scattering process simpli¬ 
fies enormously: the single quark scattering dom¬ 
inates any other subprocesses and we have essen¬ 
tially a Compton scattering on a single quark. 
Compared with DIS in which the nucleon breaks 
into many fragments, DVCS in a sense is a non- 
invasive surgery. 

The relation of the DVCS amplitude to the 
energy-momentum form factors can be seen as 
follows. The Compton scattering amplitude is a 
matrix element of the time-ordered product of the 
electromagnetic currents, 


= i 


'd^ze'^-{p'\TJ^ (-|) r {^)\p)m 


where q — {q + q')l2. In the Bjorken limit, —q^ 
and P ■ q oo and their ratio remains finite, we 
can use the operator product expansion to ap¬ 
proximate the current product, 

TJ^^iOriO) = ... + C"^"“^(e")?7a*D^V' + .-(17) 

where is the coefficient function and other 

local operators have been omitted. Therefore, the 
DVCS amplitude contains the matrix element of 
the quark energy-momentum tensor fl^ . The 
gluon contribution comes at the next-to-leading 
order. However, as the ellipses indicate, the 
DVCS amplitude contains much more informa¬ 
tion than just the energy-momentum tensor. It 
contains a whole new type of parton distributions. 


q q 




Figure 2. (a) Compton scattering on a nucleon; 
(b) Leading hand-bag contribution to DVCS. 


To motivate the distributions, let us consider a 
form factor of a local current Its matrix 

element between the nucleon states |P) and \P') 
can be obtained (at least in the infinite momen¬ 
tum frame for a good component of the current) 
by annihilating a quark of momentum k at the 
spacetime point 0 and creating the same particle 
back at the same point with a new momentum 
q + k. On the other hand, the ordinary Feynman 
parton distribution q{x) can be interpreted as a 
forward matrix element in the nucleon state in 
which a quark of momentum xp^ is annihilated 
at point 0 and created back with the same mo¬ 
mentum at a different spacetime point An^, where 
n is a light-cone vector = 0 conjugating to the 
nucleon momentum p. 

The generalized parton distributions (GPDs) 
combine the spacetime structures of the above 
matrix elements. They are defined as an inter¬ 
ference amplitude in which a quark is annihilated 
at spacetime point 0 with momentum xp^, and 
created back at another point An^ with a new 
momentum xp^ -I- q^, yielding a recoiled nucleon. 
Therefore, the GPDs naturally include the elastic 
form factors and Feynman parton distributions in 
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their kinematic limits. In fact, they contain much 
more information than the two traditional types 
of observables. 

Technically, GPDs can be defined through the 
matrix elements of a bilocal light-cone operator 

HHBl, 


1 f dX 

2 I 


e^\x / p, 


^0 i 


X Ve 


-w L7'- da n-A{an) 




P 


= -U{P') i^U{P)+Eq{x,i,t) 


2 ^ ’ 2M ^ ’ 


( 18 ) 


The light-cone bilocal operator (or light-ray op¬ 
erator) arises frequently in hard scattering pro¬ 
cesses in which partons propagate along the light- 
cone. The parton distributions are most naturally 
defined in terms of the matrix elements of the 
bilocal operator. In this context, the Feynman 
X is just the conjugating Fourier variable of the 
light-cone distance. The Lorentz structures in the 
second line in the above equation are independent 
and complete. 

The GPD’s are more complicated than the 
Feynman parton distributions because of their de¬ 
pendence on the momentum transfer A. As such, 
they contain two more scalar variables besides the 
Feynman variable x. The variable t is the usual 
t-channel invariant which is always present in a 
form factor. The ^ variable is a natural product 
of marrying the concepts of the Feynman distri¬ 
bution and form factor: The former requires the 
presence of a preferred momentum along which 
the partons are predominantly moving, and the 
latter requires a four-momentum transfer A; ^ is 
just a scalar product of the two momenta. 

Since the quark and gluon energy-momentum 
tensors are just the twist-two, spin-two, parton 
helicity-independent operators, we immediately 
have the following sum rule from the off-forward 
distributions; 


1-1 


dxx[Hq{x,^,t) + Eq{x,^,t)] 


— ^q{t) + Bq{t) , 


( 19 ) 


where the ^ dependence, or Cq{t) contamination, 
drops out. Extrapolating the sum rule to t = 0, 
the total quark (and hence quark orbital) contri¬ 
bution to the nucleon spin is obtained. A similar 
sum rule exists for gluons. Thus a deep under¬ 
standing of the spin structure of the nucleon may 
be achieved by measuring the GPDs from high 
energy experiments. 

Recently, M. Burkardt has constructed an in¬ 
terpretation of GPD in the coordinate space [|^. 
Gonsider the nucleon state localized in the trans¬ 
verse plane at r_L = 0, 

\p+,r± = 0) = N J d‘^p±\p+,p±) , ( 20 ) 

where A is a normalization factor. Define a par- 
ton distribution q(x,b±) which is the density of 
partons with longitudinal momentum xp^ and 
the transverse distance b± in the state, 

q{x,bx)=J 

X (p+,ri|V'(0,6)7+V^(An,6)|p+,r_L). (21) 

Then it can be shown that q{x, b) is the Fourier 
transformation ofH{x, — A5_, ■C = 0) with resprect 
to the transverse momentum transfer, 

q(x,b) = J d‘^AxH{x,-Al)E’’-^^ . ( 22 ) 

Thus the GPDs provide the transverse location 
of the partons in the nucleon. 

From the viewpoint of the low-energy nucleon 
structure, it is, perhaps, most interesting to con¬ 
sider GPDs as the generating functions for the 
form factors of the so-called twist-two operators. 
Recall that the matrix elements of the electro¬ 
magnetic current in the same nucleon state are 
determined by symmetry, whereas those in the 
unequal momentum states define the (Dirac and 
Pauli) form factors which contain such interesting 
information as the charge radius and magnetic 
moment of the nucleon. The following tower of 
twist-two operators represents a generalization of 
the electromagnetic current 

O"- = ^ql ■■■l , (23) 

where all indices are symmetrized and traceless 
(indicated by (• ■ •)) and V= {V — V)/2. Techni¬ 
cally, these operators transform as (n/2,n/2) of 
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the Lorentz group. They appear, as we have dis¬ 
cussed before, in the operator production expan¬ 
sion of the two electromangetic currents. Thus, 
although these generalized currents do not couple 
directly to any known fundamental interactions, 
they can nonetheless be studied indirectly in hard 
scattering processes. 

Since the operators for n > 1 are not related 
to any symmetry in the QCD lagrangian, their 
matrix elements between the equal momentum 
states, 

, (24) 


contain valuable dynamical information about 
the internal structure of the nucleon. The ^ 
dependence of the above matrix elements signi¬ 
fies the dependence on renormalization scales and 
schemes. The quark distribution q(x, intro¬ 
duced by Feynman has a simple connection to 
the above matrix elements: 


j dxx^ ^q{x,^) = an{p) 


(25) 


X P,„,2z(i)A^" ■ • ■ 

i=0 

+ C,„(t) Mod(n+ 1,2) ■Lu{P’)U{P) 

••■A'"’*) , (26) 

where U{P') and U{P) are the Dirac spinors and 
Mod(n-|-l, 2) is 1 when n even, 0 when n odd. For 
n > 1, even or odd, there are n -I- 1 form factors. 
Cqnit) is present only when n is even. 

Like the forward matrix elements a„, the above 
form factors define GPDs completely. Introduce 
a light-like vector which is conjugate to P in 
the sense that P -n = 1. Write P = p+ (M^/2)n, 
where M = — t/4 and p is another light-like 

vector. Contracting both sides of Eq. ( p6|) with 
, we have 

= U{P') fiUiP)Hq„itt) 
+UiP')^^^^^^^U{P)Eqr.{U) , (27) 

where 


where q{x, p) is chosen to have support in (—1,1). 
For a: > 0, q{x,p) is the density of quarks which 
carry the x fraction of the parent nucleon momen¬ 
tum. The density of antiquarks is customarily de¬ 
noted as q{x,p), which in the above notation is 
-q{-x,p). 

Just like the form factors of the electromag¬ 
netic current, additional information about the 
nucleon structure can be found in the form fac¬ 
tors of the twist-two operators when the matrix 
elements are taken between the states of unequal 
momenta. Using Lorentz symmetry and parity 
and time reversal invariance, one can write down 
all possible form factors of the spin-n operator 



= U{P')-f^^^U{P) 


Y 41,„.2*(<)A^= • • • ■■■P 




2 = 0 


U{P')^^U{P) 


i=0 

+Mod(n+l,2) Cqn{t){-2^T , 

[^] 

Eqni^.t) = Y 

i=0 

-Mod(n+ 1,2) C,„(t)(-2e)"(28) 

Here we have defined a new variable ^ = —n-Al2. 
Clearly, the ^ dependence of Hqn and Eqn helps 
to distinguish among the different form factors of 
the same operator. Now the moments of H{x, t) 
and E{x,^,t) are simply 



dxx'^ ^Eq{x,^,t) 

II 


dxx^~^ Hq{x,^,t) 

= Hqr,{C,t) . 

(29) 


Since all form factors are real, the new distribu¬ 
tions are consequently real. Moreover, because of 
time-reversal and hermiticity, they are even func¬ 
tions of 


2M 
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The Compton amplitude can be expressed in 
terms of the GPDs. In the leading order, 

1 1 


'jnfll' _ 


91 


dx 


X — ^ + ie X + ^ — ie 




»2 
q 




dx 


' -1 


X — ^ + ie 




»2 
9 


X + ^ — ie 

(30) 


where n and p are the conjugate light-cone vec¬ 
tors defined according to the collinear direction 
of P, and is the metric tensor in the trans¬ 
verse space. ^ is related to the Bjorken variable 
XB = —q^/{2P ■ q) by xb = 2^/(1 -I- ^). In addi¬ 
tion to DVCS, the electromagnetic radiation from 
the scattering lepton, the so-called Bethe-Heitler 
process, also produces the same final state (back¬ 
ground). 

The first evidence for DVCS was seen by HI 
and Zeus collaborations at the HERA collider, 
DESY ^ . In Fig. 3, we have shown the 
differential cross section for e^p —> e^^p as a 
function of from HI. The hatched histogram 
shows the contribution of the Bethe-Heitler pro¬ 
cess. There is a clear access of events above the 
background. 

At low-energy facilities such as JLab, the 
Bether-Heitler process overwhelms the cross sec¬ 
tion. In this case, one can isolate the interference 
contribution from the Bethe-Heitler and DVCS 
processes. Shown in Fig. 4 is the single beam 
spin asymmetry measured by CLAS as a func¬ 
tion of the angle between the virtual and outgo¬ 
ing photons 1^. A similar data exist from the 
HERMES collaboration More experiments 
on DVCS will be done in the future at HERA, 
CERN, and Jefferson Lab. 

An all-order proof of the DVCS factorization 
was discussed in and in Factoriza¬ 

tion concerns separation of soft physics at the 
scale of hadron masses and perturbative physics 
at the scale of probe in the DVCS amplitude. 
This can be done at the same level as for inclusive 
DIS. Therefore experimental cross section can be 
used to directly extract GPDs. This is nontrivial 



Figure 3. Cross section measurements for the 
1*P IP DVCS process as a function of 


beause the final photon here is on-shell. Techni¬ 
cally, it amounts to showing all soft physics either 
can be absorbed into the GPDs, which are non- 
perturbative anyway, or are down by powers of 
l/Q^. The factorization result can be formulated 
in terms of an operator product expansion famil¬ 
iar in the field theory textbooks. 

The next-to-leading order corrections to DVCS 
are important for a precision extraction of GPDs 
and for an estimate of scaling violation. The com¬ 
plete result is now known. The one-loop correc¬ 
tions to DVCS have hrst been studied by Osborne 
and Ji [^, and also by Belitsky and Mueller | p^ . 
Two-loop anomalous dimensions were obtained 
by Belitsky and Muller ||2^. It is the same as 
that for inclusive DIS except for the non-diagonal 
contribution which can be determined by the one- 
loop conformal anomaly. Classical chromody¬ 
namics is invariant under comformal symmetry 
which is broken by quantum mechanical effects. 
Because of this, an operator can acquire a scaling 
dimension under the spatial conformal transfor¬ 
mation (conformal anomaly). The size for the 
NLO corrections has been estimated in [p^ . 

The leading higher-twist contribution to DVCS 




















Figure 4. Single spin asymmetries measured by 
the CLAS collaboration which is sensitive to the 
inteference between Bethe-Heitler and DVCS pro¬ 
cesses. 


comes from twist-three longitudinal photon scat¬ 
tering. In the interference amplitude, the imagi¬ 
nary part has a characteristic sin 2(j) and real part 
cos 2^ dependence. The contribution can be es¬ 
timated in the so-called Wilczek-Wendzura ap¬ 
proximation (neglecting dynamical higher-twist 
effects). Belitsky, Mueller, and Kirchner et al. 
1^^ found that the correction to the single spin 
asymmetry from the beam polarization is about 
6%, and from the target polarization the correc¬ 
tion is about 9%. For the spin-averaged cross sec¬ 
tion it is 17%, and 3% for the double spin asym¬ 
metry. 

DVCS can be generalized to the case of exclu¬ 
sive production of mesons Collins, Frankfurt 
and Strikman have shown that the deep-exclusive 
meson production is factorizable to all orders in 
perturbation theory j^. The vector meson pro¬ 
duction entails a rich spin and flavor structure. 
For example, the vector meson production is sen¬ 
sitive to quark helicity-independent distributions, 
whereas the pseudo-scalar mesons are sensitive 
to quark helicity-dependent distributions. The 


next-to-leading order perturbative QCD correc¬ 
tion for pseudo-scalar meson production has first 
been calculated by Belitsky and Mueller |3|] . Me¬ 
son production may be easier to detect; however, 
it has a twist suppression, l/Q^. In addition, 
the theoretical cross section depends on the un¬ 
known light-cone meson wave function. Recently, 
it has been found that a quantity very sensitive to 
the quark angular momentum is the target trans¬ 
verse spin asymmetry for vector meson produc¬ 
tion p^ .The asymmetry comes from the inter¬ 
ference of two DVCS amplitudes and is linear in 
the distribution. Since the asymmetry 

involves the ratio, it is insensitive to the next-to- 
leading order and higher-twist effects. The vec¬ 
tor meson final state allows separation of different 
quark flavors. For example, the production is 
sensitive to the combination 2 -|- J^, the oj me¬ 

son is sensitive to 2J„ — Jd, and finally the 
meson is sensitive to the — Jd combination. 

To summarize, the recently discovered general¬ 
ized parton distributions are a new type of nu¬ 
cleon observables which contain rich information 
about the internal structure of the nucleon. Many 
experiments have been planned to measure the 
distributions. Lattice QCD will be a unique the¬ 
oretical tool to calculate the distribution from the 
first principles [^ . 
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